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Abstract 

Using the adiabatic switching of interactions, we estabhsh a condition for the 
existence of electronic quasiparticles in a Luttinger hquid. It involves a charac- 
teristic interaction strength proportional to the inverse square root of the system 
length. An investigation of the exact energy level separation probability distri- 
bution shows that this interaction scale also corresponds to a cross-over from the 
non interacting behaviour to a rather typical case for integrable systems, namely 
an exponential distribution. The level spacing statistics of a spin 1/2, one branch 
Luttinger model are also analyzed, as well as the level statistics of a two coupled 
chain model. 



The field of strongly correlated electron systems has recently stimulated interesting discus- 
sions which are sometimes challenging some more traditional ideas on the many body problem. 
For instance, Anderson has proposed that the low energy properties of a two dimensional Hub- 
bard model are not properly described by a Fermi liquid theory In a recent paper 0, he 
emphasizes that this question requires a non perturbative treatment, and a careful consid- 
eration of boundary conditions. As a consequence of the difficulty of the problem, a lot of 
effort has been recently dedicated to numerical investigations either with Monte Carlo meth- 
ods or exact diagonalizations However, the available sizes remain quite small, and the 
interpretation of these results is often delicate. A rather different approach has been proposed 

recently with the hope to develop new tools for extracting more information from finite 
systems studies. These works have shown that for a large class of low dimensionnal strongly 
correlated systems, the energy levels exhibit statistical properties rather well described by 
random matrix theory. For instance, a regime of energy level repulsion is clearly seen in most 
investigated cases, at the exeption of integrable models such as the nearest neighbor or the 
interaction Heisenberg spin chain. Such a behavior has been extensively discussed in the 
context of quantum chaos. More precisely, it has been verified that many time reversal sym- 
metrical classicaly chaotic systems generate a spectrum in good agreement with the Gaussian 
Orthogonal ensemble predictions |^. By contrast, simple integrable systems yield in general 
uncorrelated energy levels and the usual exponential distribution for energy level spacings [|] . 

In this paper, we are investigating some possible connections between simple physical 
properties of an interacting Fermi system, such as the existence of long lived electronic quasi- 
particles and the energy level distribution. Intuitively, if the energy levels of the interacting 
system keep a simple one to one correspondence with those of the non interacting system, we 
expect on one hand Fermi Liquid Theory to be valid, and on the other hand the statistical 
properties of the spectrum to remain qualitatively similar as for the free electron case. The 
interpretation of random matrix behaviour is not straightforward. It may simply indicate that 
a model is non integrable. For a normal Fermi Liquid, electronic quasiparticle are expected 
only at low energies compared to the Fermi energy. Furthermore, already for the particle-hole 
phase spaces, interactions induce new collective modes, such as the zero sound, and the idea of 
a one to one correspondence with the non interacting gas does not hold for the whole spectrum. 
Clearly, it would be very interesting to see if the spectrum of a normal Fermi Liquid exhibits 
some features which would distinguish it from a random matrix Hamiltonian. However, this 
would likely require an intensive numerical effort (since best candidates would be at least 
two dimensionnal systems). For the sake of simplicity, and the motivation of doing analytical 
calculations, we have concentrated in this work on a one dimensional model, namely the Lut- 
tinger model [0 , which is integrable at any coupling strength. Interestingly, this feature holds 
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for any system length [^]. Furthermore, it provides a good example of a non- Fermi liquid, 
which can be viewed as a non translation invariant fixed point for many interacting systems 
in one dimension. 

This paper is organized as follows. A first part investigates the condition for the exis- 
tence of electronic quasiparticles, using the adiabatic generation of eigenstates. An existence 
condition is established, from the combined requirement of having a negligeable generation of 
non adiabatic components and absence of decay. This criterion is satisfied if the interaction 
strength is less than a constant divided by the square root of the system length. As expected, 
no quasiparticles are found for an infinite system at any finite value of the interaction parame- 
ter. This result is also rederived from a simple a nalysis of the single particle Green's function 
for a finite system. The second part is devoted to the study of the level spacing distribution 
as the interaction is gradually increased. We show that the typical interaction scale locating 
the departure from the highly degenerate non interacting system towards a more generic inte- 
grable model with a Poisson distribution is the same as the previous one. So, for this simple 
situation, noticeable change in the energy level distribution is reflected by the disappearance 
of electronic quasiparticles. Then, the last two sections of this paper are dedicated to variants 
of this model, namely in the spin 1/2 case and forward scattering only, for both one and two 
coupled one dimensionnal systems. A brief conclusion summarizes our results. 



1 Adiabatic switching on of interactions 

A formal way to generate quasiparticles in an interacting Fermi liquid is to apply the Landau 
switching on of interaction procedure, namely to start from a free particle added above the 
Fermi sea, and to switch on interactions adiabatically. The corresponding time dependent 
hamiltonian is: 

H = Ho + Voe'\ (1) 

where the interactions term Vq is switched with a rate e. 

Provided it is successfuU, this procedure establishes a one to one correspondance between 
the free gas excitations, and the dressed excitations of the Fermi hquid, namely, the quasipar- 
ticles. For a Fermi liquid, the validity condition of this procedure is 0: 

r(efc) < e < efc, (2) 

where is the energy of the quasiparticule, with respect to the Fermi surface, and r(efc) is the 
decay rate of the quasiparticule. For a normal Fermi liquid, one can show p that r(efc) ~ e|. 
At small energies, r(efc) ^ e^, so that it is possible to choose a rate e to perform the switching 
on procedure. 
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The aim of this section is to investigate under which conditions the switching on procedure 
is vahd in a one dimensional Luttinger hquid. We shall henceforth exhibit an inequality similar 
to equation (|^) for the rate e in the case of a Luttinger liquid. 



1.1 Introduction 

We first wish to sum up some results concerning the formalism developped in |^. This will 
also permit us to fix the notations, which shall be used in the rest of the paper. 

The fermions are on a ring of perimeter L, with periodic boundary conditions, so that the 
wave vectors are quantized {k = ^n, with n an integer). 

As we treat only low energy properties of a spinless, one dimensional Fermi gas, the cur- 
vature of the dispersion relation may be neglected. The two linear branches in the dispersion 
relation emerging from each extremity of the Fermi surface are extended to arbitrary energies. 
This linearized model is the Luttinger gas model, which hamiltonian is: 

= Vf ^{pk - kp) : 4pCfcp :, (3) 

kp 

where vp is the Fermi velocity and p = +1 or —1 labels the branch (right or left). We shall 
also use the real space field ippix) associated to the right (left) free fermions. Furthermore, 
c^p is the Fourier transform of ^pp{x)■. 

Notice that the sign of the phase factor is not arbitrary, but is chosen such as right moving 
fermions with a positive wave vector propagate to the right. 

Because of the presence of an infinite number of fermions in the ground state, the density 
operators 

PqP ~ 5Z • ^k+q,p^k,p '■ (5) 
k 

have anomalous commutation relations (Schwinger terms): 

[Pgp, P-q'p'] = — ^^pp'^gq'- (6) 

They may consequently be used to build a set of boson creators (? 7^ 0). To handle the 
real space bosonic field, one needs to define 

%{x)=p^N,-tY^eipq)i^y/'e^^^a,. (7) 

^ g^O Mil 
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The q = modes correspond to charge and current excitations. Their algebra involves the 
unitary ladder operators Up constructed in p . They act only in the q = sector, and increase 
by one the charge on the p branch. The complete form of the bosonic fields, including the 
q = modes, is: 

^p(x) =^p + $p(x) + <l'+(x), (8) 

where 9p is the phase conjugate to Np. 

We shall also use the important relation to pass from a real space boson description to a 
real space fermion description: 



Expressed on this new basis, the free hamiltonian becomes: 

H^ = vfY1 I^Ks + VFjiNl + Nl), (10) 

where Nji (N^) denote the number of right (left) moving fermions added above the vacuum 
state. In terms of charge N = + A^^^ and current J = — Nl variables, the energy of the 
charge and current excitations is: Vf^{N'^ + J^). 

Note that the action of the boson creation operators and of the ladder operators on the 
ground state generates a basis of the Hilbert space. The completeness may be shown by 
comparing the generating functions of the degeneracies (ie the finite temperature partitions 
functions) for both the free electrons basis and the boson basis. The notation for the kets of 
the second basis is: 

We now briefly describe the formalism to deal with interactions. The two-particle interac- 
tions term is written as: 

= 7 E Kp.pP-.-p- (12) 
^ pq 

For simplicity, our treatment does not include interactions between fermions lying on the 
same side of the Fermi surface. Only g2 interactions are relevant in the physics we shall 
develop. One important feature of the interactions Vq is that they are cut off for impulsions 
greater than the inverse of a length scale R. We shall use the following expression of Vq (for 
q < l/R): 

Vq = V{l-{qRr). (13) 
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The intensity of the interactions is parametrized by V, and the shape of Vq is parametrized 
by a. The bosonized form of the interaction Hamiltonian is: 



ZL IL 

The total hamiltonian is diagonalized by the following Bogoliubov transformation: 

6^ = coshi^qO^ — sinh(y9qa_q, (15) 
where the angle ipq is defined as: 

Va 

tanh2(/?g = ^. (16) 

Vp 

The total hamiltonian reads, after the Bogoliubov transformation: 

The effect of the interactions is to give a non zero ground state energy: 

Eo = Wiujq-VFq), (18) 
^ q 

where 

^q = {vl-v;f'w. (19) 

Interactions also shift the energies of the oscillators from VF\q\ to ujq. Finally, charge and 
current excitations acquire different velocities vn = vse~'^'^ and vj = vse'^'^. In these relations. 
If is the infrared limit of ipq and the sound velocity vs is related to the infrared limit of the 
dispersion relation (|T9|): 

vs = Umq_^oivl~V^'Y/'. (20) 

In the presence of interactions, one needs to normal order the field ipp{x) in terms of 6^ 
bosons, which leads to the following expression of ^p{x) : 

Tlx 

<^p{x) = p—Np -i^{9{pq) cosh - e{-pq)sinh<fq)e"^''bq. (21) 
The fermion field reads, in terms of bosons: 

<?>o 
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1.2 Interaction picture for c^ji\{Np}) 

As t — > — oo, the system is made up of a right moving fermion, with an impulsion k added above 
a Dirac sea \{Np}), and interactions are vanishing. This section deals with the propagation of 
this state, clj^\{Np}) , as interactions are switched on. 

The first step is to decompose the state c^fi\{Np}) into bosonic modes. The action of c'^j^ 
on the vacuum |{iVp}) in the q = sector is simply to increase by one the number of right 
moving fermions, by the action of the ladder operator Ur. 

To obtain the action of c^^ in the q ^ sectors, we first Fourier transform c^^ into the 
real space field ^^(a;) for right-moving fermions. 

We replace the expression of 0p (x) in (|^) by its expression (|^) in terms of bosonic modes 
a^. The developpement of the exponential e"**^*^^^ leads then to an expression of cl^\{Np}) 
as a linear combination of bosonic states, with occupation numbers {rig}: 



4r\{Np}) = j:4j2<l^<i-ik-{kF + ji2Nn+m] (23) 
n^fyr) ' \{Nn + l,N,},{ng}). (24) 

q>0 \/nql K^QJ 

The delta function insures that only bosonic states with a total impulsion equal to A; — 
kp — j^{2Nf> + 1) survive in the decomposition. As no interaction couples the two branches, 
creating a right-moving fermion does not generate left moving bosons. 

The second step is to propagate the bosonic wave packet (pS]) . Instead of dealing with the 
rather complicated superposition ( p3D of bosonic states, we focus on the propagation of a single 
term IjA'p}, {uq}). We shall use the bosonized form of the two-body interaction hamiltonian, 
and look for a solution of the time dependent Schrodinger equation: 

. ^i^M^'<^K HL\{N^,MUW. (25) 
The "int" label stands for an interaction picture. The initial conditions are: 

limt-.-oo\{Np}, K}),„t(t) = \{Np}, {ug}). (26) 

The bosonic states are propagated under the form of a coherent state: 

\{Np}, {ng})Ut) = iV(^,(t))e-^«"^>'*) n e-'^^^'^^>-'>\{Np}, {nj). (27) 

q>0 
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The prefactor N{z) normalizes |{A^p}, {nq})int{t): 



= n (l - k.P) ' ■ (28) 

To determine the time dependent (j){{nq},t) and {zq(t)} functions, we first change Zq{t) into 
Uq(t), with Zq(t) = 'Ug(t)e^™^'^*, and then identify both sides of the Schrodinger equation. We 
obtain first order non hnear differential equations for {uq{t)} and (j){{nq},t): 



"^""'^^^ '-2tVFqUq{t) = qVqm-u'Jt)) (29) 



dt 

d(j){{nq},t) 



J2inq + l)qVqit)Imiuqit)). (30) 



We have discarded in equation |3^ a term depending only on and J, which leads only 
to a global phase factor. Translated in terms of and Zq variables, the initial conditions 
(|26|) simply mean that (j){t) and Zq{t) are vanishing as t ^ — oo. These differential equations 
describe the propagation of a single component of the wave packet (^). The propagation of 
the summation is obtained as a superposition of the different components after propagation: 



{4R\{N,})Uit) = Y.s{j:Qriq-ik-{kF + ji2Nn + m] (31) 

M \i>o ^ I 

riq 

\{Nn + hN,},{nq})Ut). (32) 



q>0 \ friq\ 



1.3 Adiabaticity condition 

We are looking for a solution of equation (|SD| ) which depends only on the variable s = et, in the 
small e limit. It is possible since the external time dependance in equation ( |30D involves only 
et. We assume then Uq{s) = n°(s) + €u^{s) + O(e^). Neglecting the O(e^) terms in equation 
(|5DD leads to: 



2tVFqu'^^is) = qVqis){l-uf\sf) (33) 

(m) -)- 2 

ds 



dvP 

-^{s) + 2iVFqu^^\s) = -2qVq{s)u^^\s)u^^\s), (34) 



where: 



Vq{s) = \/V. (35) 
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The purely adiabatic solution u'^\s) is given by: 



Kis) = t/z-A-vf + ^Jvl - V,{sf) = ^tanhy,;(.). 



(36) 



Using this solution in equation (0) gives the first finite e correction: 




0, which 



(37) 



(38) 



2q{vl - V^) 



This condition depends explicitely on g, and is satisfied for any value of q if 



e -C 4:Ttvf/L. 



(39) 



Here, we assume a weak coupling, namely \Vq\ ^ vp- It should be noticed that this 
upper bound on e is a much more restrictive condition than the corresponding upper bound 



quasiparticules of the Landau theory are not exact eigenstates of the interacting system. 
They are obtained in a situation where the thermodynamic limit is taken first, whereas the 
generation of exact eigenstates would require e to go to zero as the typical spacing between 
energy levels. Our criterion (|39| ) corresponds to this second situation. This choice has been 
motivated by the possibility to construct the exact eigenstates of a Luttinger liquid. 

1.4 Adiabatic propagation in a Bogoliubov subspace 

The aim of this section is to propagate a fermion during the switching on procedure. We 
suppose that the condition (|3^) is satisfied, and we now look for a minoration of e. We first 
search an approximation for the evolution operator in the limit e <^ x^^' order 
e°, the evolution operator ?7e(0, — oo) realizes the Bogohubov transformations of angles {^g}, 
corresponding to the rotation of the basis of eigenstates as interactions were switched on from 
zero at time t = —oo to {v^g} at time t = 0. We shall note f/°(0, —oo) the corresponding part 
of the evolution operator. [7° must have the property that: 




We interpret this as a consequence of the fact that the 




(40) 



This equality is verified if U has the following form: 




(41) 



q>0 
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To see it, we differentiate each operator U^a^{U^) ^ and U^a^q{U^) ^ with respect to ^9° and 
solve the differential system. 

However, at higher orders in e, the evolution operator must take into account the phase 
factor (/)({?T,g}, t), the evolution of which is given by the equation (|30|). Assuming that the 
propagation is adiabatic, we approximate ImUq{t) in ( pO| ) by ImvP^it): 

Imuq{t) ~ tanhv9j(s = et). (42) 

We use the expression (|^) for V^(t), and integrate the differential equation (|30|) for the 
phase factor 0({ng},t). A constant (infinite) phase factor associated to the propagation of 
the ground state is factored out. Thus, we obtain the form of the evolution operator in the 
adiabatic limit (at order e° for the operator ?7°, and at order 1/e for the phases): 



U^-oo) = t/°exp.{E^(v^°)n (43) 

q>0 ^ 

= exp{E^;(a:-a^^-Sa_,)}exp{z5:^(^;)n. (44) 

q>0 q>0 

In the integrations, we have assumed that the interactions are weak, and the phase factors are 
given, at the lowest order in 99°. 

The rest of this section is devoted to the calculation and the interpretation of the overlap: 

F{x-x')r^F{x,x') = {{Np}\^R{x')Ur\0,-oo)^U^)U,i^,-oo)\{Np}), (45) 

between the dressed fermions \l/^(x)f/e(0, — cx))|{iVp}), and the bare ones: \E'^(x')f/~^(0, — oo)|{A'p}). 
To perform it, we use the expression ( p2D of the field for right moving fermions, and the ap- 
proximation ( |4^ for the evolution operator. The computation is straightforward, and F{x, x') 
is the product of three terms: 

1) a phase term 

N = exp {-i{kp + ^{2Nr + l))(a; - x')}, (46) 

corresponding to the propagation in the g = sector. 

2) A term corresponding to the left moving bosons normal ordering in (^): 

Gi = exp{-5:|^(sinh^;)2} (47) 

q>0 ^1 

3) A term coming from the right moving bosons normal ordering: 

G2(x, x') = exp ^e-*^(^-^'-^^^)} (48) 
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The result for the overlap is: 



F{x, x') = yNGiGiix, x'). (49) 



The Gi term contains the usual physics of the orthogonality catastrophy [jlOl . If (/?° is assumed 
to be constant between q = ^ and q = l/R, and zero afterwards, and if L 3> -R, Gi can be 
calculated as: 

^ ^|)-sinhV° (50) 

In the weak coupling limit, one can deduce the characteristic interaction scale associated to 
the orthogonality catastrophy: 

^o.c. = VF{\n^)-'/'. (51) 

To obtain the energy scale associated to the G2 term, we use the relation (|13D and approximate 
the phase as: 

e Avpe 2vFt 

The first term is linear in q up to the impulsion scale l/R. The second term is associated to 
smaller impulsion scales. The formers are relevant for a quasiparticle. If k is the impulsion 
of the quasiparticle with respect to the Fermi level, the energy scale Vdeph associated to the 
dephasing is given by: 

-^{kRT = 27r, (53) 

ZVp€ 

that is: 

The switching on procedure shall henceforth be successfuU provided the intensity of interac- 
tions V is much smaller than Vdephik), that is: 

kV^{kRY 

4:71V I 

1.5 Conclusions 



< e. (55) 



For the switching on procedure to create a quasiparticle, the conditions ( P^ and ( ^3]) have to 
be simultaneously satisfied, that is: 

kV^{kRY ATI . ^, 

— — < e < —vp. (56) 

AiTVp L 

This inequality is satisfied if the following consistency condition is fuUfilled: 

Ativp 



10 



As we shall see, this condition has a simple interpretation on the spectrum of the Luttinger 
model. At this stage, we should again emphasize that the upper bound on e is more restrictive 
than in Landau theory. If we use the more usual condition that the spread in energy is smaller 
than the average density of the wave packet, equation (p6D is replaced by: 

kV\kRY 



Attvj 

and the consistency condition is: 



< e < kvF, (58) 



V « ( j^Y'W. (59, 

The absence of Landau quasiparticule in the thermodynamic limit is then atributed to orthog- 
onality catastrophy, as indicated by equation (^). 

1.6 Comparison with the Green's function 

In this section, we calculate the Green's function for the finite size Luttinger model: 



GR{x,t;x',t') = -z{({iV,}|e*^(*'-*Vi?(a;')e^'''^*'~*VA(a;)|{^^p})^(^'-i) (60) 
-(x ^x';t^ t')} 

and reestablish the consistency condition (|57|). Note that in equation (|BDD, \{Np} > denotes 
an eigenstate of the interacting system. 

To calculate the Green's function, we use the expression (^ of the field ipni^) terms of 
the Bose field, and normal order the expression ( |60D of the Green's function with respect to 
the bosonic modes 6^. The computation is straightforward, and the result is: 

GR{x,t;x',t') = Z!e^(fcF+-A)(-'-)e^f(-iv(27v+i)+.,(2j+i))(t'-t) (g^^ 
exp (-2^— (sinh ^^f) 
{ [exp (^ ( ^ ) (cosh (^,) 2e^^("' -^'^ e"'"' (*' "*) ) 

g>0 

exp (^(— ) (sinh (^,)2e-'''("'-")e-'"'(*'-*))] 

-[x^x']t^t']e{t-t')} 

The dispersion in the frequencies leads to decoherence after a time t^. {k is the impulsion 
of the quasiparticle, with respect to the Fermi level), tk may be estimated in the same way as 
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we did for Udeph, and one finds: 

= V^kiTur (62) 

For a system of size L, the wave packet is stable, provided it can cross the ring without 
decoherence: 



vptk > L, (63) 

that is: 

Up to some numerical dimensionless constants, this criterium is the same as the consistency 
condition (0) for the switching on of interactions. 



2 Level statistics of the interacting Luttinger model 
2.1 Introduction 

We first need to find out a proper sector of the Hilbert space, in which we shall compute the 
level statistics. We note Hj^ the subspace with given current J and charge N. 

In the free case, the boson basis of Hjn can be organized as follows: consider all the sets of 
occupation numbers {n^} such as, for all q, rig = or n° ^ = 0. The corresponding states 
are annihilated by any pair destruction operator: aga_g|{n°}) = 0. Starting from and 
creating pairs generates a subspace Hpairs{{ng}) . A basis of impairs ({^°}) is made up of all the 
states |{n° +p\q\}) with arbitrary occupation numbers for the pairs {pq}q>o- H^j is the direct 
sum of all the impairs ({'^g}) subspaces. 

The subspaces Hpairs{{n^q\) remain stable under the action of the interaction hamiltonian 
H^, so that they are appropriate to the study of the levels evolution. 

We choose N = J = Q and drop the energy term associated to {^g}, since we always handle 
differences between consecutives levels. The energy levels are given by: 

E(K}) = 5:2w,(l-(^)^)i/2 (65) 
where we use the expression (|13D for Vq. 
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2.2 Description of the algorithms 

In this section and the next paragraph, we use reduce units for the energies and impulsions: 
cu is an energy divided by and q is an impulsion divided by 
The degeneracies of the Luttinger model are given by: 

= (66) 



L/2 fj^ 1 

^ n -^e— ^ (67) 



Replacing the S function by its integral representation leads to: 

Let g^^^Lj) be the number of different sets of occupation numbers, having the property that: 

oJ^f^lni. (68) 

l=k 

Of course, g'^^^co) — g{oj). The integral representation for g^^\u) reads: 



Using the integral representations for g^^\uj)^ we obtain the following recurences: 

^W(u;) = X:^^^n^-^), (70) 

i'=k 

which allows us to numerically compute g{u!). 

With a similar recursion, we may generate all the states of the free Luttinger model: the 
states with an energy cu are obtained by adding a boson with an impulsion u on the states 
with an energy cu — u. 

As far as the interacting Luttinger model is concerned, we need to generate all the energy 
levels with an energy inferior as a given cut-off cuo- Since there are an infinite number of 
levels in the sector under consideration, we need to introduce such a cut-off to compute the 
statistics. We shall then compute the statistical properties of this set of levels. If a sufficient 
number of levels with an energy inferior as ujq has been generated, the statistical properties are 
independant on cuq. To generate the levels, we remark that the frequencies of the oscillators 
increase with their impulsion. So that we successively fill up the individual oscillator levels, 
starting with the smallest frequencies. 
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2.3 Level statistics 

2.3.1 Degeneracies of the free Luttinger model 



Using the recursion relation ([70|) , we computed the degeneracies of the first 800 levels of the 
free Luttinger model. The asymptotic form of the density of states may be derived in terms of 
initial fermions. The partial degeneracies for n-particules n-holes excitations in a one branch 
model are: 

9^''\^)= E E S{u-j:u{h)-±u{k'^). (71) 

{ki}i^l...„ {k'.}i^i,,,„ 1=1 i=l 

The sets {ki} ({A;^}) are the impulsions of the holes (particules) , and are constrained by the 
Pauli principle ^ kj {k[ ^ k'-) for all indices i ^ j- This sum is approximated by assuming a 
constant density of states, neglecting the Pauli exclusion principle, replacing the discrete sum 
by an integral: 

g^^^Uj) = -—r / duJi / duJ2--- / duJ2nS{uJ - {uJi + ... +UJ2n)) (72) 



(n!)2 Jo Jo Jo 
The multiple integral is readily evaluated and leads to: 

+00 +00 2n—l 

.M^E.<">(.)^E(;,5,p^. (73) 

For sufficiently large energies, the sum may be approximated by its saddle point value, ap- 
proximately reached for the following value of n: 

^* = (74) 

The degeneracy evaluated at n = n* is: 

23/4 ^ 

~ (27r)3/2 ^5/4 ^^P^8u;- (75) 

We computed the summation ([73| ) in order to test the accuracy of the saddle point approxi- 
mation, which is plotted on figure |1|. The exact degeneracies of the Luttinger model reveal to 
be inferior as the saddle point asymptotic form, which is imputed to the exclusion principle 
(figure p. 

2.3.2 Qualitative structure of the spectrum 

The evolution of some energy levels as a function of the interactions is plotted on figure ^. In 
this spectrum, we distinguish two regions: 
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1) No level crossings are present at sufficiently small energies and interactions. The free 
Luttinger model (V = 0) belongs to this part of the spectrum. In this region, the statistics 
are ill defined for they strongly depend on the energy cut offi 

2) If E and V are large enough, level crossings occur, and level statistics are Poisson 
statistics. The convergence of the statistics as a function of the energy cut off cq is shown 
on figure |^. Here, we emphasize that these level crossings occur because the Luttinger model 
remains integrable at any value of the coupling constant. 

To characterize the separation between these two regions of the spectrum, the location of 
the crossings is estimated in the following way: as the intensity of interactions V is equal to 
zero, the spectrum is made up of equidistant degenerate levels, separated by an amount of 
energy AE = ^vp- As V is turned on, the degeneracies are lifted. We focus on a single fan 
of levels. All the levels are degenerate if V = 0, and their energy is E^ = 2vFk, where k is 
the total impulsion of the states. For a given value of V, all the levels lie between Emm and 
Emax- Emin is obtained as all the quanta are in the smallest energy state (namely q = ^vp), 
so that: 



Emin = 2vFk{l - (76) 

Vp 

Emax corresponds to a state with one quantum in the highest q = k state: 

Emax = 2VFk{l-^f'\ (77) 

As the interaction parameter V increases, the levels evolve and the first crossings occur as the 
width of the fan Emax — Emm is of order AE. This condition defines the interaction energy 
beyond which crossings exist: 

Via bosonization, the free Luttinger liquid is described as a set of harmonic oscillators 
with commensurable frequencies. As interactions are switched on, the oscillator frequencies 
vary and become incommensurable. In |0, Berry and Tabor show that a system with a 
finite number of generic harmonic oscillators does not exhibit level clustering. It appears that 
increasing the number of oscillators with incommensurable frequencies generates clustering. 

2.3.3 Quasi particle destruction and level spacing statistics 



The condition ( |75D separates two regions of the spectrum. The same energy scale controls the 
existence or the absence of a quasiparticule in a Luttinger liquid, in the sense of adiabatic 
continuation of exact eigenstates. We have thus shown that the structure of the spectrum of 
the finite size Luttinger liquid is related to the succes or the failure of adiabatic generation of 
eigenstates from the non interacting fermion system. 



15 



2.3.4 Limit R = 



Consider the case of the two branch Luttinger hquid with = +00, and a constant interaction, 
namely, for all q, Vq = V. In that case, all the bosonic modes keep their coherence whatever 
the value of e. The decoherence time tk, given in (|62D, is infinite. The condition ( ^5]) associated 



to the dephasings is always verified whatever the value of e. The only remaining restriction 
for the switching on procedure to be successful is thus: 

.«^. (79) 

The level statistics are singular in this limit. The degeneracies of the fan of levels are never 
lifted, whatever the intensity of interactions V. The degenerate levels depend on V in the 
following way: 

^(K}) = E2w.(l-(-)')^/'. (80) 

However, we note that the overlapp between the eigenstate thus constructed and the state 
obtained from the action of the bare electron operator on the interacting ground state is 
vanishing according to equation (|^) since R = 

3 Level spacing statistics for a spin 1/2, one branch Lut- 
tinger model. 

The rest of the article is devoted to the study of some models derived from the two branch, 
spinless Luttinger liquid model. We begin with the one branch Luttinger model, with spin 
1/2, and a g4, interaction. The kinetic energy term is: 



H'' = Y,VF{k-kF):ctck„:, (81) 

fccr 

where the label a denotes the spin component along the z axis. The interaction is given by: 



The usual spin and charge combinations: 

(P.T+P.i) (83) 

{Pq]-Pqi). (84) 

have bosonic commutation relations, and the total hamiltonian H = + is diagonal in 
terms of spin and charge variables: 





(- 


J 1/2 








s: = 


(f 


J 1/2 




Lq 
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H = vcY: qC^C, + vsY: qS^S, + vp^{N^ + Nf). (85) 

q>0 q>0 '^^ 

The charge and spin velocities are: fc = "^f + |^ and vs = vp — 
The g4 interaction is switched on adiabatically: 

94{t) = gV. (86) 

The evolution operator is: 

f/e(0, -oo) = exp {-% ^ iz^q{ncq - ^Sg)}, (87) 

where ncq = C^Cg and nsq = S^Sg. The overlap 

F{x,x') = {{Np}\^^{x')U~\0,-oo)^:^{x)U,{0,-oo)\{Np}) (88) 

is found to be equal to: 

F{x, x') = ^e^(f (2^T+i)+fc^)(--') 1 . 1 . . (89) 

^ (1 _ e^¥(^-^'-&))i/2 (1 _ e*¥(^-^'+&))i/2 

Spin charge separation is effective if the real space separation is of order which leads to 
the energy scale for spin charge decoupling: 

9l = ^^ (90) 

where k is the impulsion of the quasiparticle with respect to the Fermi surface. The switching 
on procedure is sucessful provided <^ gl- Since the transformation (|^) is independent on 
the interactions, there is no upper limit for the rate of switching on e. 

In the same way as for the Luttinger liquid, the sector of the Hilbert space has to remain 
stable under the action of the evolution operator (|^). Since t/e(0, — oo) is diagonal in term of 
charge and spin variables, the relevant sector has a given impulsion k. This sector corresponds 
to a single fan of levels, with no crossings, except for (74 = 0, leading to singular statistics. One 
may compute the statistics in the whole Hilbert space, namely to superpose the uncorrelated 
blocs with different impulsions. The statistics still remain singular. The degeneracies of some 
levels are not lifted for any value of the interaction (74. These singularities correspond to re- 
maining degeneracies as the impulsions of the charge and spin part are specified independently, 
and are reminiscent of the degeneracies of the free Luttinger model. The spectrum exhibits 
further singularities at non zero level spacings, due to the linear dependence of the levels in 
5^4: the statistics do not become poissonian even though uncorrelated sectors are superposed. 
Note that many degeneracies, and the singularities at non zero level spacings are expected to 
disappear if (74 is not a constant as a function of q. In this more generic case, the Poisson 
statistics is expected. 
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4 Level spacing statistics for a model of 2 coupled chains 



We now discuss the level statistics for a model of two coupled Luttinger liquids. This model 



is solved in [|TT| and we first remind some results. 
The two chains kinetic energy is given by: 

^° = ^FE(^-M:c+aA. :, (91) 

kcxa 

where a labels the chain and a the spin. The interactions consist of a (74 term: 

= § E P^-P^.-., (92) 

kaa 

and of a hopping term between the two chains: 

H^ = -t^T.4a.Ck-aa. (93) 
kcra 

Only the case of two coupled one branch models is treated. This is sufficient since no inter- 



action couples right and left fermions. Fabrizio and Parola [11| were able to diagonalize the 



hamiltonian H = H^ + H'^ + H-^. The excitation spectrum of the model exhibits four branches: 

ep(g) = Upq (94) 
e^(g) = u^q (95) 

e+(g) = liup + Ua)q + J {^{up - u^)qy + 4tl (96) 



= l{'^P + Ua)q-]J{\{up-u^)qy + Atl (97) 

The ground state is such as all the states with a negative energy are occupied, and all 
the states with a positive energy are empty. We computed the level statistics for a toy model 
with only the €-{q) branch, in a sector of given total impulsion q. We study the evolution of 
the statistics as the dimensionaless hopping constant t_L = is fixed, and = varies. 
The statistics exhibit a cross-over between two regimes as decreases. This cross-over is 
controlled by the same lenght scale ^ = "^""'^ as in |ll[]. If g,^ <^ 1, the statistics are singular, 



with a sharp peak at s = 0. In this regime, the dispersion relation e_(g) may be approximated 
as: 

e-iq) = ^iup + u^)q-2t^. (98) 

The linear q dependance induces high degeneracies in the excitation spectrum, leading to a 
sharp peak for zero separation. 
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In the opposite regime {q^ ^1), the statistics are poissonian. The corresponding spectrum 
is plotted on fig. ^ In that case, the curvature of the dispersion relation e_(g) is no longer 
negligeable, and individual fermion levels can no longer be considered as equidistant. 

Notice that the cross-over observed here is similar to the case of the one dimensional, one 
branch Luttinger liquid with g-dependant interactions. In both cases, the dispersion relation 
is linear as the interaction parameter is set to zero (corresponding to a highly degenerate spec- 
trum), and becomes non linear as interactions are switched on (leading to a random spectrum). 
This transition is independent of the bosonic or fermionic nature of the particules. In the one 
dimensional Luttinger liquid, we dealt with bosons, and the particules under consideration in 
the case of the two coupled chains are fermionic. 

What happens if we now take the four branches into account? In the regime ^ 1, we 
observe a peak for s = 0, coexisting with a poissonian distribution for non zero separations (see 
fig. 1^, where the peak is suppressed for clarity). The peak for s = is due to the degeneracies 
in the excitation spectrum, induced by the presence of the two linear branches. An exemple 
of such degenerate configurations, with 2 particule hole excitations is as follows: the two holes 
have impulsions hi and /i2, and belong to the e_ branch. The particules with impulsions pi 
and p2 are on the linear ep branch. Consider an other excitation, deduced from the previous 
one as follows: the holes have the same impulsions {h[ = hi and /ig = h2). The impulsions of 
the particules are such as pi+ P2 = p'l + p'2- Since all the particules belong to the same linear 
branch, these configurations are degenerate. 

Thus, the existence of the two regimes q^ ^ 1 and ^ 1 in the coupled chains is re- 
fiected in the statistical properties of the spectrum. To summarize, we have studied a special 
class of models, since they are integrable for any value of the coupling constant. In general, a 
non interacting fermionic quasiparticle can be described as a linear combination of degenerate 
eigenstates, which undergo an energy splitting as interactions are switched on. This is respon- 
sible for the decay of such a quasiparticle state, and provides a lower bound for the switching 
rate e, in the process of adiabatic construction of quasiparticles. The same degeneracy lifting 
has been found to modify the energy level spacing distribution, from a singular behaviour for 
a degenerate, non interacting system, to a more generic Poisson distribution already observed 
in many integrable systems. We should stress that both aspects are non universal features of 
the models. More precisely, they depends on the complete g-dependance of the interaction 
functions g2 and g^. By contrast, universal properties such as correlation function exponents 
depend only on the q = limit of the couplings. We have seen that the vanishing of the 
quasiparticule residue, due to orthogonality catastrophy is also such an universal property, 
independant on the fine structure of the spectrum and its statistics. 

In this paper, we couldn't address the question of strongly correlated fermion systems 
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leading to gaussian orthogonal ensemble (G.O.E.) statistics. However, the present work indi- 
cates that one of the most interesting questions is whether the difference between G.O.E. or 
Poisson distribution is a universal feature of a low-energy fixed point or not. Our paper has 
been dedicated to fine tuning phenomenas within an integrable class of models, and the lack 
of universality found here is not surprising. Intuitively, the difference between Poisson and 
G.O.E. statistics is much more robust and might still be a way to distinguish between several 
physically non equivalent fixed points. 

R. M. whishes to thank J.C. Angles d'Auriac for help with programmation, P. Degiovanni 
for discussions about bosonisation and acknoledges the hospitality of NEC Research Institute 
where part of this job was performed. 



20 



References 

[1] P.W. Anderson Phys. Rev. Lett. 64, 1839 (1990); Phys. Rev. Lett. 65, 2306 (1990) and 
Phys. Rev. Lett. 66, 3226 (1991). 

[2] P.W. Anderson Phys. Rev. Lett. 71, 1220 (1993). 

[3] For a recent review, see E. Dagotto, National High Magnetic Field Laboratory preprint 
(1993), Florida State University 

[4] G. Montambaux, D. Poilblanc, J. Behssard and C. Sire, Phys. Rev. Lett. 70, 497 (1993); 
T.C. Hsu and J.C. Angles d'Auriac, Phys. Rev. B 47, 14291 (1993); D. Poilblanc, T. 
Ziman, J. Belissard, F. Mila and G. Montambaux, preprint (1993). 

[5] A review on the connections between random matrix theory and chaos is by O. Bohigas, 
in "Chaos and Quantum Physics", proceeding of the Les Houches Summer School of 
Theoretical Physics, Les Houches, 1989, edited by M.J. Giannoni, A. Voros and J. Zinn- 
Justin (North Holland, Amsterdam, 1991). 

[6] M.V. Berry and Tabor, Proc. R. Soc. Lond. A 356, 375 (1977). 

[7] J.M. Luttinger, J. Math. Phys., 4, 1154 (1963). 

[8] F.D.M. Haldane, J. Phys. C, 14, 2585 (1981). 

[9] P. Nozieres, "Le probleme a N corps", Dunod, Paris 1963. 

[10] P.W. Anderson, Phys. Rev. Lett. 18, 1049 (1967). 

[11] M. Fabrizio and A. Parola, Phys. Rev. Lett 70, 226 (1993). 



21 



Figure captions 



Figure 1: 

Degeneracies of tlie free Luttinger model, compared to tlie saddle point approximation. 
{log5f(ci;) — I log 2 + |log27r + ^\og^/uJ}/^/8uJ is plotted as a function of ^/uJ. This func- 



tion equals 1 for the saddle point approximation. In plot (1), g{uj) is the exact degeneracies. 
As expected, the saddle point approximation overevaluates the degeneracies since it takes into 
account particule-hole exitations forbidden by the exclusion principle. In plot (2), all the 
terms of the summmation ([73|) are taken into account. The saddle point approximation in 
([73|) underevaluates the degeneracies, and becomes exact at high energies. 



Figure 2: 

Evolution of some levels as a function of interactions, is a linear decreasing function of 
q, such as (pq=o = 2.5a, and (Pq>26L/2n = 0. a parametrizes the interaction strenght, and the 
energy is in units of vp^. For the plot to be readable, all the levels are not shown. 



Figure 3: 

Evolution of the level spacing statistic as a function of the cut-off Cq. fq is a decreasing linear 
function, such as 95g=o = 0.25 and (p26^ — 0- The statistics converge slowly to a Poissonian 
distribution ("expo"). The statistics are plotted for cq equal to 5, 9, 13. The number of levels 
taken into account in the statistics is repectively: 4196, 97438, 1048214. 



Figure 4: 

Level spacing statistics for the model of two coupled chains in the regime ^ 1. Only 
excitations of the lowest energy branch e_(g) are taken into account, and the analysis is 
restricted to the Ip — Ih and 2p — 2h excitations only, for parameters equal to: p = 200, 
tper = 20, g = 0.5, with the following notations: p is the total impulsion divided by 
iper = and g = ir^- 24000 states were generated. The value of the parameter is 100. 



Figure 5: 

Level spacing statistics for the model of two coupled chains in the regime 3> 1, with the 
four excitation branches. Only \p — Ih and 2p — 2h excitations were taken into account. 
The parameters are set to: p = 100, tper = 100, g = 0.5 and qC, = 50. The number of 
computed levels is 42692. Among them, 8293 separations are equal to zero. For visibility, the 
level statistics is cut off for separations inferior as 0.05, which supresses the large pic at zero 
separations. 
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